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Abstract 

Arguments are summarized, that neutral matter made of helium, carbon, 
etc., should form a quantum liquid at the above-atomic but below- nuclear 
densities for which the charged spin-0 nuclei can condense. The resulting sub- 
stance has distinctive features, such as a mass gap in the bosonic sector and 
a gap-less spectrum of quasifermions, which determine its thermodynamic 
properties. I discuss an effective field theory description of this substance, 
and as an example, consider its application to calculation of a static potential 
between heavy charged impurities. The potential exhibits a long-range oscil- 
latory behavior in which both the fermionic and bosonic low-energy degree of 
freedom contribute. Observational consequences of the condensate for cooling 
of helium-core white dwarf stars are briefly discussed. 
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Foreword 



Like many in the audience, I first met Misha on the pages of journal pubhcations, 
before meeting him in person. While working on an undergraduate thesis at Moscow 
University, I came across Misha' s review paper "Anomalies and Low- Energy Theo- 
rems of Quantum Chromodynamics" [1]. Impressions of that work were very distinct 
- a clear exposition of subtle field theory aspects of the quantum anomalies, cul- 
minating in creative applications to low-energy hadron phenomenology. The work 
stood out by its originality, depth, inspiration and balance of the formalism and 
applications - the remarkable signatures of Misha's enormous contribution to theo- 
retical physics at the forefront of both field theory and particle phenomenology. 

I met Misha in person in Minneapolis in 1998. The discussion with him was very 
inspiring. Soon, in Aspen, we started to work on a project. A bit later I ceased the 
opportunity to get exposed to two years of a unique FTPI experience. We continued 
to work on and off on various projects since then. I value those works very highly, 
and feel privileged, as I'm sure many of you do too, for having such a collaborator. 

Happy 60th Birthday Misha! 
Description of charged condensate 

Consider a neutral system of a large number of nuclei each having charge Z, and 
neutralizing electrons. If average inter-particle separations in this system are much 
smaller than the atomic scale, ~ 10~^ cm, while being much larger than the nuclear 
scale, ~ 10~^'^ cm, neither the atomic nor nuclear effects will play any significant 
role. Moreover, the nuclei can also be treated as point-like particles. 

In what follows we focus on spin-0 nuclei with Z < 8 (helium, carbon, oxygen), 
and consider the electron number-density in the interval Jo ~ (0.1 — 5 MeVY. Then 
the electron Fermi energy will exceed the electron-electron and electron-nucleus 
Coulomb interaction energy. Moreover, at temperatures below ~ 10"^ K, which are 
of interest here, the system of electron represents a degenerate Fermi gas. 

Since the nuclei (we also call them ions below) are heavier, temperature at which 
they'll start to exhibit quantum properties will be lower. Let us define the "critical" 
temperature Tc, at which the de Broghe wavelengths of the ions begin to overlap 



where, ttih denotes the mass of the ion (the subscript "H" stands for heavy) , and d 
denotes the average separation between the ions^. 

"'^The de Broglic wavelength above is defined as XdB = 27r/|k|, where k^/2m^f = 3kBT/2. We 
define Tc as the temperature at which A^b — d. Note that this differs by a numerical factor 
of i/27r/3 from the standard definition of the thermal de Broglie wavelength, A = ^/2n/mkBT, 
that appears in the partition function of an ideal gas of number-density n in the dimensionless 
combination A^n. 
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Somewhat below quantum-mechanical uncertainties in the ion positions be- 
come greater than an average inter-ion separation. Hence the latter concept looses 
its meaning as a microscopic characteristic of the system; the ions enter a quantum- 
mechanical regime of indistinguishability. Then, the many-body wavefunction of the 
spin-0 ions should be symmetrized, and this would unavoidably lead to probabilistic 
"attraction" of the bosons to condense, i.e., to occupy one and the same quantum 
state. We refer the system of condensed nuclei and electrons as charged condensate. 

In the condensate the scalars occupy a quantum state with zero momentum. 
Moreover, small fluctuations of the bosonic sector happen to have a mass gap, 

= {Ze^jQ/mnY^'^i which exceeds by more than an order of magnitude. There- 
fore, once bosons are in the charged condensate, their phonons cannot be thermally 
excited. However, the gap-less fermionic degrees of freedom are thermally excited, 
and carry the most of the entropy of the entire system [2]-[5]. 

For further discussions it is useful to rewrite the expression for Tc in terms of the 
mass density p = tuhJo measured in g/cm^: 

where the baryon number of an ion was assumed to equal twice the number of 
protons, A = 2Z (true for helium, carbon, oxygen...). Thus, for p = 10^ g/cm^ and 
helium-4 nuclei we get Tc ~ 10^ K, while for the carbon nuclei with the same mass 
density ~ 2 • 10^ K. 

Temperature at which the condensation phase transition takes place, Tcondens, 
need not coincide with Tp. Moreover, we would expect Tcondens ^ ^c- Calculation 
of Tcondens from the fundamental principles of this theory is hard. However, we can 
obtain an interval in which Tcondens should fit. For this we introduce the following 
parametrization: 

Tcondens — C -^c ) (3) 

where C is an unknown dimensionless parameter that should depend on density more 
mildly than does. Numerically, however, this parameter should vary in the in- 
terval 0.1 ^ C ^ 1- The point C = 0.1 would corresponds to the temperature of 
the Bose- Einstein (BE) condensation of a free gas for which, T^o^^ens — l-3/mij(i^, 
is known from the fundamental principles. The condensation temperature in our 
system should be higher than T^^^^ since the repulsion makes easier for the con- 
densation to take place [6]. In our case, repulsive interactions between the bosons 
are strong - the Coulomb energy is at least an order of magnitude greater that any 
other energy scale in the system. Hence, we should expect ( ^ 0.1. On the other 
hand, given the definition of the parameter C cannot be greater than unity. In 
what follows we will retain ( in our expressions, but use C — 1 when it comes to 
numerical estimates. 

The condensation will take place after gradual cooling, only if Tcondens is greater 
than the temperature at which the substance could crystallize. A classical plasma 



3 



crystallizes when the Coulomb energy becomes about ~ 180 times greater than the 
average thermal energy per particle [7, 8, 9]. This gives the following crystallization 
temperature^ 



Note that the density dependence of Tc is different from that of Tcryst - for higher 
densities Tc grows faster, making condensation more and more favorable! One can 
define the "equality" density for which Tcondens — T'cryst^ 



For helium, Z — 2, and peq — 10^ g/cw?, while for carbon, Z = 6, and peq — 
10^ gjcrr? (as mentioned above, we use C, ~ 1). These results are very sensitive to 
the value of for instance, pcq could be an order of magnitude higher if C — 0.5. 
Irrespective of this uncertainty, however, the obtained densities are in the right 
ballpark of average densities present in helium-core white dwarfs ~ 10^ 9 /cm, (for 
carbon dwarfs, they're closer to those expected in high density regions only [5].) 

Is the charged condensate a ground state of the system at hand? For the higher 
values of the density interval considered, the crystal would not exist due to strong 
zero-point oscillations. At lower densities, the crystalline state has lower free energy 
(at least near zero temperature) due to more favorable Coulomb binding. Hence, 
the condensate can only be a metastable state. The question arises whether after 
condensation at ~ Tcondens the system could transition at lower temperatures ~ Tcryst 
to the crystal state, as soon as the latter becomes available. 

In the condensate, the boson positions are entirely uncertain while their mo- 
menta equal to zero. In order for such a system to crystallize later on, each of the 
bosons should acquire energy of the zero-point oscillations of crystal ions. As long 
as this energy, ~ (Ze^ Jo/mi?)^/^, is much greater than Tcryst, no thermal fluctua- 
tions can excite the condensed bosons to transition to the crystalline state. The 
latter condition is well-satisfied for all the densities considered in this work. There 
could, however, exist a spontaneous transition of a region of size Rc to the crystal- 
lized state via tunneling. The value of Rc, and the rate of this transition, will be 
determined, among other things, by tension of the interface between the condensate 
and crystal state, which is hard to evaluate. However, for estimates the following 
qualitative arguments should suffice: the height of the barrier for each particle is 
(Ze^ Jo/mj^)^/^ = m^, while the number of bosons in the Rc region ~ R^Jq/Z. 
Hence, the transition rate should scale as exp{—m^JoRc/Z). Since we expect that 
Rc > 1/m-y, the rate is strongly suppressed for the parameters at hand. 

^The presented formula for the crystallization temperature is entirely classical. The temperature 
scale that determines the classical versus quantum nature of the crystallization transition is the 
Debye temperature 0£) ~ 4 • lO'^p^/^ K. Often, 6d niay significantly exceed Tcryst [10]. In such 
cases, quantum zero-point oscillations should be taken into account. This seems to delay the 
formation of quantum crystal, lowering Tcryst from its classical value at most by about ~ 10% [11]. 
Since this is a small change, we will ignore it in our estimates. 




(4) 




(5) 



4 



Effective field theory description 

We use a low-energy effective field theory description to study the charged conden- 
sate. Even though realistic temperatures in the system may be well above zero, 
we focus on the zero-temperature limit. The relevance of this limit is justified 
a posteriori and goes as follows: the spin-0 nuclei undergo the condensation to 
the zero-momentum state; their phonons cannot be excited since their gap, m^, is 
greater than T^. On the other hand, gap-less near-the-Fermi-surfacc quasielectrons 
will be excited. Therefore, all the thermal fluctuations will end up being stored in 
the fermionic quasiparticles. For the latter, however, the finite temperature effects 
aren't significant since their Fermi energy is so much higher, T/jJ^ < IQ-^. We 
note that the finite temperature effects, in a general setup with condensed bosons, 
were calculated in Refs. [12, 13]. 

We begin at scales that are well below the heavy mass scale m/^, but somewhat 
above the scale set by max[/x/,me], where fif and rUe are the electron chemical 
potential and mass respectively. Hence the electrons are described by their Dirac 
Lagrangian, while for the description of the nuclei we will use a charged scalar 
order parameter As it was shown in [4], in a non-rclativistic approximation 

for the nuclei, the effective Lagrangian proposed by Greiter, Wilczek and Witten 
(GWW) [14] in a context of superconductivity, is also applicable here, given that an 
appropriate reinterpretation of its variables and parameters is made. 

The construction of the GWW Lagrangian is based on the following fundamen- 
tal principles: it is consistent with the translational, rotational, Galilean and the 
global U{1) symmetries, preserves the algebraic relation between the charged current 
density and momentum density, gives the Schrodinger equation for the order param- 
eter in the lowest order, and is gauge invariant [14]. Combined with the electron 
dynamics the GWW Lagrangian reads (we omit for simplicity the Maxwell term): 

A// = V (^|(rL>o$ - {D,^m - + i^ii^Dl - ruf)^, (6) 

where we use the standard notations for covariant derivatives with the appropriate 
charge assignments: Dq = {do — iZeAo), Dj = (dj — iZeAj), Dj[ = (9^ + ieA^^, while 
V{x) stands for a general polynomial function of its argument. The coefficients of 
this polynomial, V{x) = ^n^^i ^-re dimensionful parameters that are inversely 

proportional to powers of a short-distance cutoff of the effective field theory^. 

^In general one should also add to the Lagrangian terms unr^*^, A($*$)^/m|^, 

Xi[^*^)iptp/[mHJo ), and other higher dimensional operators that arc consistent with all the 
symmetries and conditions that lead to (6) (the Yukawa term is not). Here fx^R denotes a 
non-relativistic chemical potential for the scalars. These terms are not important for the low- 
temperature spectrum of small perturbations we're interested in. as long as A, Ai ^1 and Jq <C m^. 
However, near the phase transition point it is the sign of /Ltjvfi that would distinguish between the 
broken and symmetric phases, so these terms should be included for the discussion of the symmetry 
restoration. We also note that the scalar part of (6) is somewhat similar to the Ginzburg-Landau 
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Once the basic Lagrangian is fixed, we introduce the electron chemical potential 
term /ifip^ip , where /if = ep = [(Stt^ Jq)^/"^ + m^]^/^. This is the only term that at 
the tree level sets a frame in which the electron total momentum is zero. 

There exists a homogeneous solution of the equations of motion that follow from 
the effective Lagrangian (6) [3]: 

Z|<l>|2 = ./o, A^ = 0, V'iO) = l. (7) 

(We use the unitary gauge $ = |$|). The condition V'{0) = 1 is satisfied by any 
polynomial functions V{x) for which the first coefficient is normalized to unity 

V{x)^x + C2x'' + .... (8) 

The above solution describes a neutral system of negatively charged electrons of 
charge density — eJo, and positively charged scalar condensate of charge density 
Ze$+$ = eJo [4, 5]. 

Calculation of the spcctrTim of small perturbations is straightforward. The La- 
grangian density for the fluctuations in the quadratic approximation reads [2] 

A = -^F^l + - \m-^A] + \ ^i^I^^Ao , (9) 

where A denotes the Laplacian, and the last term emerged due to mixing of Aq with 
the fluctuation of the |$|, which we integrated out. As before, 

m = , (10) 

^ rriH 

and ttIq = m^, + 6*26^ Jq. At this stage we retained the fermionic fluctuations only in 
the Thomas- Fermi approximation [3]; an important reflnement of this approxima- 
tion, discussed in [4], will be included below. 

That there are no pathologies in (9), such as ghost and/or tachyons, can be seen 
by calculating the Hamiltonian density: 

^ = f + ^ + + {d^^^) + \m',A]. (11) 

Here, = run'm^ and ttj = —Fqj. The Hamiltonian is positive semi-definite. 
Moreover, the spectrum has a mass gap determined by (10). There are two 
transverse polarizations of a massive photon, as well as the longitudinal mode, the 
phonon, with the same mass [2]. 



(GL) Lagrangian for superconductivity. However, there are significant differences between them, 
one such difference being that the coherence length in the GL theory is many orders of magnitude 
greater than the average interelectron separation, while in the present case, the "size of the scalar" 
$ is smaller that the average interparticle distance. 
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The massive bosonic collective excitations give rise to exponentially suppressed 
contributions to the value of specific heat of the charged condensate since typically 
ruj ^ Tc- The suppression scales as exp(— m^/T), where T ^T^. This is in contrast 
with the crystal, where the dominant contribution to the specific heat comes from 
a gap-less phonon, and scales with temperature as T^. 

As to the electrons, their behavior is similar in both crystal and condensate 
cases. At temperatures of interest they form a degenerate Fermi gas with gap- 
less excitations near the Fermi surface. Their contribution to the specific heat scales 
linearly with temperature. In the case of crystallized substance this is sub-dominant 
to the specific heat due to the crystal phonon. For the charged condensate, however, 
the (quasi) electron fiuctuations are the dominant contributors to the specific heat. 

To study the cfi'ccts of collective bosonic and fermionic modes, as an interesting 
example, we look at a potential between two impurity nuclei (say hydrogen, or 
helium-3) of charge Qi and Q2. The calculation of the propagator that involves the 
hght collective modes (for relativistic fermions) gives the following result [4] : 

^ ^ fe'^"^ , 4aemfcFsin(2fcFr)\ , , 

Vstat = «e.Qig2 (^^cos(Mr) + J ^ ^s/ ) " (12) 

The first, exponentially suppressed term modulated by a periodic function, is due 
to cancellation between the screened Coulomb potential and that of a phonon [4]. 
The fact of such a cancellation, and that it could give rise to the oscillatory behavior 
of the exponentially screened potential was pointed out before in Ref. [15] in the 
context of superconductivity^. 

Most important, however, is the second term in (12) that has a long-range [4]. 
It dominates over the exponentially suppressed term in (12) for scales of physical 
interest, and exhibits the power-like behavior modulated by a periodic function. 

The potential (12) is a generalization of the Friedel potential to the case when 
in addition to the fermionic excitations there are also collective modes due to the 
charged condensate. The long-range oscillating term in (12) is also a result of a 
subtraction between the conventional Friedel term and the long-range oscillating 
term due to a phonon. As a result, its magnitude is suppressed compared to what 
it would have been in a theory without the condensed charged bosons [4] (see, [16] 
for the discussion of the conventional Friedel potential, and Ref. [13] for its recent 
detailed study in the presence of the charged condensate at finite temperature.)^ 

The potential (12) is not sign-definite. In particular, it can give rise to attrac- 
tion between like charges; this attraction is due to collective excitations of both 
fermionic and bosonic degrees of freedom. This represents a generalization of the 
Kohn-Luttinger [18] effect to the case where on top of the fermionic excitations the 
collective modes of the charged condensate are also contributing^. 

""■I'd like to thank Ki-Myeong Lee who recently brought the paper [15] to my attention. 
®Note that for spin-dependent interactions the same effects of the charged condensate would 
give a generalization of the Ruderman-Kittel-Kasuya-Yosida (RKKY) potential [17]. 

^In the charged condensate Cooper pairs of electrons can also be formed, however, the cor- 
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Applications to White Dwarfs 

The above described system of electrons and nuclei constitutes cores of white dwarf 
stars. Up to a factor of a few, these arc roughly Earth's size solar-mass objects; 
their mass density may range over ~ (10® — 10^°) g/cm?, most of them being near 
the lower edge of this interval. Since the dwarf stars exhausted thermonuclear fuel 
in their cores already, they evolve by coohng [19]; the ones that we consider in this 
work cool from ~ 10^ K down to lower temperatures. 

As a typical dwarf star cools down, the Coulomb interaction energy in a classical 
plasma of charged nuclei will significantly exceed their classical thermal energy, and 
the nuclei, in order to minimize energy, would organize themselves into a crystal 
lattice [20] . In most of these cases quantum effects of the nuclei should be negligible; 
for instance, the Debye temperature should be less than the temperature at which 
crystallization takes place, and the de Broglie wavelengths of the nuclei should be 
much smaller than the average internuclear separations. This indeed is the case in 
majority of white dwarf stars, the cores of which are composed of carbon and/or 
oxygen nuclei and span the interval of mass densities around ~ (10^ — 10^) g/crn^. 

However, there exists a class of dwarf stars in which the nuclei enter the quantum 
regime before the classical crystallization process sets in [10, 11]. Among these, 
furthermore, there is a relatively small subclass of the dwarf stars for which the 
temperature Tc, is higher than the would-be crystallization temperature Tcryst [5]. 
In such dwarf cores the charged condensation should be expected to take place. 

White dwarfs composed of helium constitute a smaller sub-class of dwarf stars 
(see, [21, 22] are references therein); they exhibit best conditions for the charged 
condensation. Most of helium dwarfs are believed to be formed in binary systems, 
where the removal of the envelope off the dwarf progenitor red giant by its binary 
companion happened before helium ignition, producing a remnant that evolves to a 
white dwarf with a helium core. Helium dwarf masses range from ~ 0.5 Mq down to 
as low as (0.18 — 0.19) Mq, while their envelopes are mainly composed of hydrogen. 

Using the approach of [23], and following [5] we will consider an over-simplified 
model of a reference helium star of mass M = 0.5 Mq with the atmospheric mass 

responding transition temperature, and the magnitude of the gap, are suppressed by a factor 
exp(— l/Cgjj), where e^^j is proportional to the value of the inter-electron potential that contains 
both screened Coulomb and phonon exchange. The fact that this potential has attractive domain, 
but is very small, is suggested by the static potential found in [4] (see also eq. (12) above); the 
latter is suppressed by a power of a large scale M. Furthermore, taking into account the frequency 
dependence of the propagator in the Eliashbcrg equation docs not seem to change qualitatively the 
conclusion on a strong suppression of the Green's function and pairing temperature. 

Hence, even though the bosonic sector (condensed nuclei) is superconducting at reasonably high 
temperatures <^ 10® K, interactions with gap- less fcrmions could dissipate the superconducting 
currents. Only at extremely low temperatures, exponentially close to the absolute zero, the elec- 
trons could also form a gap leading to superconductivity of the whole system. In the present work 
we consider temperatures at which electrons are not condensed into Cooper pairs, and ignore the 
finite temperature effects. 
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fractions of the hydrogen, and heavy elements (metalhcity) respectively equal to 

X ~ 0.99, Zm ^ (0.0002 - 0.002) . (13) 

The lower value of the metalhcity ~ 0.0002 is appropriate for the recently 
discovered 24 He WDs in NGC 6397 [22] , but for completeness, we consider a wider 
range for this parameter. 

It is straightforward to find the following expression for the cooling time of a 
star in the classical regime [23] 



tne ^ 



CAm 



u 



3-5 5 TT^/Td-S 3' 

-{TP - TP) + Z-^^{Tp - TP) 



(14) 



where Tf and Tq denote the final and initial core temperatures. The first term 
in the bracket on the right hand side corresponds to cooling due to classical gas 
of the ions and the second term corresponds to the contribution coming from the 
Fermi sea. The latter is sub-dominant in the range of final temperatures we are 
interested in (the factor Z in front of this term is due to Z electrons per ion) . Since 
Tf <^ To, the age of a dwarf star typically doesn't depend on the initial temperature. 
Neglecting the fermion contribution, we find time that is needed to cool down to 
critical temperature Tf = T^ 

Where an order of magnitude interval in (15) is due to the interval in the envelope 
metalhcity composition given in (13). We also find the corresponding luminosities 

M fT \ '^^^ 

L(T,) ^ (10« ergis)— (^^J ^ 1.5 • (10"^ - 10-^)Le , (16) 

which are in the range of observable luminosities {Lq ~ 3.84 • 10^^ erg/s). 

After the condensation, specific heat of the system dramatically drops as the 
collective excitations of the condensed nuclei become massive and "get extinct". A 
contribution from the Fermi sea, which is strongly suppressed by the value of Fermi 
energy, becomes the dominant one. The phase transition itself would take some time 
to complete, and the drop-off in specific heat will not be instantaneous. 

In the zeroth approximation, we can regard the transition to be very fast, and 
retain only the fermion contribution to specific heat below Tc- Then, the expression 
for the age of the star for Tf < Tc, reads as follows 



CAm 



u 



3 TT^ kn -ii -3 

-(Tc ' -Tn ') + Z- -(T, ^ - To n 

5 ^ ^ ' 3 eP ^ ° ^ 



(17) 



Notice the difference of (17) from (14) - in the former Tf < Tc and it is Tf that 
enters as final temperature in the fcrmionic part, while Tc should be taken as the 
final temperature in the bosonic part. 
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From the ratio of ages, rj = tHe/t'He^ identical helium dwarf stars, with 

and without the interior condensation, we deduce that the charged condensation 
substantially increases the rate of cooling- the age could be twenty times less than 
it would have been without the condensation phase [5]. 

The condensation of the core would induce significant deviations from the classi- 
cal curve for helium white dwarfs. What is independent of the uncertainties involved 
in these discussions, is the fact that the luminosity function (LF) will experience 
a significant drop-off after the charged condensation phase transition is complete. 
This is due to the "extinction" of the bosonic quasiparticles below the phase transi- 
tion point. In fact, the LF will drop by a factor of ~ 200. This may be relevant for 
an explanation of the observed termination of a sequence of the 24 He WD's found 
in [22]. See Ref. [5] for more details. 

Finally, the magnetic properties of the charged condensate, which are similar to 
those of type II superconductor, and in particular admit the presence of Abrikosov's 
vortices, were studied in Ref. [24]. As was shown there, only very strong magnetic 
fields, > 10^ Gauss, will be able to penetrate the dwarf cores in the vortices, while 
weaker fields will be entirely expelled from it. 
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